Large-scale in-silico analysis of CSF dynamics within the subarachnoid space of the optic nerve.

BACKGROUND: Impaired cerebrospinal fluid (CSF) dynamics is involved in the pathophysiology of neurodegenerative diseases of the central nervous system and the optic nerve (ON), including Alzheimer's and Parkinson's disease, as well as frontotemporal dementia. The smallness and intricate architecture of the optic nerve subarachnoid space (ONSAS) hamper accurate measurements of CSF dynamics in this space, and effects of geometrical changes due to pathophysiological processes remain unclear. The aim of this study is to investigate CSF dynamics and its response to structural alterations of the ONSAS, from first principles, with supercomputers. METHODS: Large-scale in-silico investigations were performed by means of computational fluid dynamics (CFD) analysis. High-order direct numerical simulations (DNS) have been carried out on ONSAS geometry at a resolution of 1.625 µm/pixel. Morphological changes on the ONSAS microstructure have been examined in relation to CSF pressure gradient (CSFPG) and wall strain rate, a quantitative proxy for mass transfer of solutes. RESULTS: A physiological flow speed of 0.5 mm/s is achieved by imposing a hydrostatic pressure gradient of 0.37-0.67 Pa/mm across the ONSAS structure. At constant volumetric rate, the relationship between pressure gradient and CSF-accessible volume is well captured by an exponential curve. The ONSAS microstructure exhibits superior mass transfer compared to other geometrical shapes considered. An ONSAS featuring no microstructure displays a threefold smaller surface area, and a 17-fold decrease in mass transfer rate. Moreover, ONSAS trabeculae seem key players in mass transfer. CONCLUSIONS: The present analysis suggests that a pressure drop of 0.1-0.2 mmHg over 4 cm is sufficient to steadily drive CSF through the entire subarachnoid space. Despite low hydraulic resistance, great heterogeneity in flow speeds puts certain areas of the ONSAS at risk of stagnation. Alterations of the ONSAS architecture aimed at mimicking pathological conditions highlight direct relationships between CSF volume and drainage capability. Compared to the morphological manipulations considered herein, the original ONSAS architecture seems optimized towards providing maximum mass transfer across a wide range of pressure gradients and volumetric rates, with emphasis on trabecular structures. This might shed light on pathophysiological processes leading to damage associated with insufficient CSF flow in patients with optic nerve compartment syndrome.

Terabyte-scale supervised 3D training and benchmarking dataset of the mouse kidney.

The performance of machine learning algorithms, when used for segmenting 3D biomedical images, does not reach the level expected based on results achieved with 2D photos. This may be explained by the comparative lack of high-volume, high-quality training datasets, which require state-of-the-art imaging facilities, domain experts for annotation and large computational and personal resources. The HR-Kidney dataset presented in this work bridges this gap by providing 1.7 TB of artefact-corrected synchrotron radiation-based X-ray phase-contrast microtomography images of whole mouse kidneys and validated segmentations of 33 729 glomeruli, which corresponds to a one to two orders of magnitude increase over currently available biomedical datasets. The image sets also contain the underlying raw data, threshold- and morphology-based semi-automatic segmentations of renal vasculature and uriniferous tubules, as well as true 3D manual annotations. We therewith provide a broad basis for the scientific community to build upon and expand in the fields of image processing, data augmentation and machine learning, in particular unsupervised and semi-supervised learning investigations, as well as transfer learning and generative adversarial networks.

Large-scale morphometry of the subarachnoid space of the optic nerve.

BACKGROUND: The meninges, formed by dura, arachnoid and pia mater, cover the central nervous system and provide important barrier functions. Located between arachnoid and pia mater, the cerebrospinal fluid (CSF)-filled subarachnoid space (SAS) features a variety of trabeculae, septae and pillars. Like the arachnoid and the pia mater, these structures are covered with leptomeningeal or meningothelial cells (MECs) that form a barrier between CSF and the parenchyma of the optic nerve (ON). MECs contribute to the CSF proteome through extensive protein secretion. In vitro, they were shown to phagocytose potentially toxic proteins, such as α-synuclein and amyloid beta, as well as apoptotic cell bodies. They therefore may contribute to CSF homeostasis in the SAS as a functional exchange surface. Determining the total area of the SAS covered by these cells that are in direct contact with CSF is thus important for estimating their potential contribution to CSF homeostasis. METHODS: Using synchrotron radiation-based micro-computed tomography (SRµCT), two 0.75 mm-thick sections of a human optic nerve were acquired at a resolution of 0.325 µm/pixel, producing images of multiple terabytes capturing the geometrical details of the CSF space. Special-purpose supercomputing techniques were employed to obtain a pixel-accurate morphometric description of the trabeculae and estimate internal volume and surface area of the ON SAS. RESULTS: In the bulbar segment, the ON SAS microstructure is shown to amplify the MECs surface area up to 4.85-fold compared to an "empty" ON SAS, while just occupying 35% of the volume. In the intraorbital segment, the microstructure occupies 35% of the volume and amplifies the ON SAS area 3.24-fold. CONCLUSIONS: We provided for the first time an estimation of the interface area between CSF and MECs. This area is of importance for estimating a potential contribution of MECs on CSF homeostasis.

High-Throughput Lossy-to-Lossless 3D Image Compression.

The rapid increase in medical and biomedical image acquisition rates has opened up new avenues for image analysis, but has also introduced formidable challenges. This is evident, for example, in selective plane illumination microscopy where acquisition rates of about 1-4 GB/s sustained over several days have redefined the scale of I/O bandwidth required by image analysis tools. Although the effective bandwidth could, principally, be increased by lossy-to-lossless data compression, this is of limited value in practice due to the high computational demand of current schemes such as JPEG2000 that reach compression throughput of one order of magnitude below that of image acquisition. Here we present a novel lossy-to-lossless data compression scheme with a compression throughput well above 4 GB/s and compression rates and rate-distortion curves competitive with those achieved by JPEG2000 and JP3D.

Synchronisation through learning for two self-propelled swimmers.

The coordinated motion by multiple swimmers is a fundamental component in fish schooling. The flow field induced by the motion of each self-propelled swimmer implies non-linear hydrodynamic interactions among the members of a group. How do swimmers compensate for such hydrodynamic interactions in coordinated patterns? We provide an answer to this riddle though simulations of two, self-propelled, fish-like bodies that employ a learning algorithm to synchronise their swimming patterns. We distinguish between learned motion patterns and the commonly used a-priori specified movements, that are imposed on the swimmers without feedback from their hydrodynamic interactions. First, we demonstrate that two rigid bodies executing pre-specified motions, with an alternating leader and follower, can result in substantial drag-reduction and intermittent thrust generation. In turn, we study two self-propelled swimmers arranged in a leader-follower configuration, with a-priori specified body-deformations. These two self-propelled swimmers do not sustain their tandem configuration. The follower experiences either an increase or decrease in swimming speed, depending on the initial conditions, while the swimming of the leader remains largely unaffected. This indicates that a-priori specified patterns are not sufficient to sustain synchronised swimming. We then examine a tandem of swimmers where the leader has a steady gait and the follower learns to synchronize its motion, to overcome the forces induced by the leader's vortex wake. The follower employs reinforcement learning to adapt its swimming-kinematics so as to minimize its lateral deviations from the leader's path. Swimming in such a sustained synchronised tandem yields up to [Formula: see text] reduction in energy expenditure for the follower, in addition to a [Formula: see text] increase in its swimming-efficiency. The present results show that two self-propelled swimmers can be synchronised by adapting their motion patterns to compensate for flow-structure interactions. Moreover, swimmers can exploit the vortical structures of their flow field so that synchronised swimming is energetically beneficial.

Mesh-particle interpolations on graphics processing units and multicore central processing units.

Particle-mesh interpolations are fundamental operations for particle-in-cell codes, as implemented in vortex methods, plasma dynamics and electrostatics simulations. In these simulations, the mesh is used to solve the field equations and the gradients of the fields are used in order to advance the particles. The time integration of particle trajectories is performed through an extensive resampling of the flow field at the particle locations. The computational performance of this resampling turns out to be limited by the memory bandwidth of the underlying computer architecture. We investigate how mesh-particle interpolation can be efficiently performed on graphics processing units (GPUs) and multicore central processing units (CPUs), and we present two implementation techniques. The single-precision results for the multicore CPU implementation show an acceleration of 45-70×, depending on system size, and an acceleration of 85-155× for the GPU implementation over an efficient single-threaded C++ implementation. In double precision, we observe a performance improvement of 30-40× for the multicore CPU implementation and 20-45× for the GPU implementation. With respect to the 16-threaded standard C++ implementation, the present CPU technique leads to a performance increase of roughly 2.8-3.7× in single precision and 1.7-2.4× in double precision, whereas the GPU technique leads to an improvement of 9× in single precision and 2.2-2.8× in double precision.